AD-A214 


O 


GL-TR«@9-C1 22 


MODTRAN: 

A  MODERATE  RESOLCTXOBS  MODEL  FOR  ifOWTRAHI  7 


Alexander  Berk, 

Lawcancfi  S.  Bernstein,  And  David  C.  Robertson 
Spectral  Sciences,  Inc, 

99  South  Bedford  Street *  #7 
Burlington,  MA  ©1803 


30  April  1989 


Final  Report 

12  May  1985  -  11  April  1989 


Approved  for  public  release; 


die tributt on 


unlim 


"This  technical  report  hat>  reviewed 


Contract  Manager 

MiKosphoric  Effects  Branch 

Optical / Inf r *r»d  Technology  Division 


and  is  approved  for  publication*’ 


Atmospheric  Effect®  Branch 

Optica  1/Inf rasred  'technology  Division 


FOR  TRK  COMMANDS® 


Gptical/Infrared  Technology  Division 


This  report  ha®  been  reviewed  by  the  ESD  Public  Affairs  Qffifla  (*A>  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) . 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
technical  information  Center.  All  others  should  apply  to  the  National 
Technical  Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  AFGL/DAA,  Hanscam  AFB,  MR  01731.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or 
notices  on  a  specific  document  require#  that  it  be  returned. 


UNCLASSIFIED 


rtTusrf  v  mwiWTiSFRjrrHirwsBr 


j 


REPORT  DOCUMENTATION  PAGE 


,  la.  REPORT  SECURITY  CLAS$IFICATi6n 
j  UNCLASSIFIED 

lb.  RESTRICTIVE  MARKINGS 

2*  SECURITY  CLASSIFICATION  AUTHORITY 

J  .  DISTRIBUTION 
Approved 

/AVAILABILITY  OF  REPORT 

"or  public  release; 

.on  unlimited 

j  7b  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

d is  tr ibut 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

j  SSI-TR-154 

S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

GL-TR-89-0122 

Ga  NAME  OF  PERFORMING  ORGANIZATION 

Spectral  Sciences,  Inc. 

6b.  OFFICE  SYMBOL 
(It  applicable) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

Geophysics  Laboratory/OPE 

1  6c.  ADDRESS  (Cl y.  Stale,  end  ZIP  Code) 

99  South  Bedford  Street,  #7 
|  Burlington,  MA  01803-5169 

7b.  ADDRESS  (City,  Start,  and  ZIP  Code) 

Air  Force  Systems  Command 

Hanscom  AFB ,  MA  01731-5000 

'  Be.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

1 

8b.  OFFICE  SYMBOL 
(It  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F19628-86-C-0079 

J  3c.  ADDRESS  (City,  Stele,  end  ZIP  Code)  j 

_ 1 

10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO. 

62101F 

PROJECT 

NO 

7670 

TASK 

NO. 

09 

U  TITLE  (Induct*  Security  Ckmfkttlon) 


MODTRAN:  A  Moderate  Resolution  Model  for  LOWTRAN  7 


1 7  PERSONAL  AUTHOR(S) 

Alexander  Berk,  Lawrence  S.  Bernstein,  and  David  C.  Robertson 


j  13a  TYPE  OF  REPORT 

Il3b.  TIME  COVERED 

14  DATE  OF  REPORT  (Year,  Month,  Dey) 

IS.  PAGE  COUNT 

|  Final 

|  FROM86MAY12  T089APR11 

89APR30 

42 

16  SUPPLEMENTARY  NOTATION 


w 

COSATI  COOES 

IB.  SUBjECT  TERMS  (Continue  on  reverse  If  neceotry  end  identify  by  block  number ) 

|  FIELD 

GROUP 

SUB-GROUP 

Transmission,  Radiation,  Infrared,  Band  Model, 

Atmospheric,  Transmittance,  and  LOWTRAN 

_ 

19  ABSTRACT  [Continue  on  reverie  H  rtecesury  end  identify  by  block  number) 


This  final  report  describee  a  new  band  model  computer  code  which  calculates 
atmospheric  transmittance  and  background  radiance  and  has  a  apectral  resolution  from 
50  to  2  cm~*.  This  new  code  is  called  MODTRAN.  The  LOWTRAN  transmittance  model  and 
associated  subroutines  have  beers  fully  integrated  into  MODTRAN,  and  all  the  LOWTRAN 
7  options,  i.e.,  aerosol  models,  atmospheric  paths,  transmitcance/radiance, 
user-specified  data,  etc.,  have  been  maintained.  Band  model  parameters  have  been 
formulated  from  the  HITRAN  line  atlas  for  twelve  atmospheric  molecules: 


h2o,  co2,  o3,  n2o,  CO,  ch4,  o2,  no,  so2,  no2,  nh3  &  hno3* 
They  were  calculated  for  i  cm  1  bins  from  0-17,900  cm-*  and  at  five  temperatures 
from  200  to  300°K.  These  parameters  are  stored  on  an  external  data  file  (tape) 
Which  is  accessed  by  the  program.  The  transmittance  is  calculated  with  an 
equivalent-width  formulation  that  takes  into  account  the  finite  spectral  width  of 
each  interval  and  the  finite  number  of  lines  contained  therein.  The  Curtia-Godson 


70  DISTRIBUTION  1  AVAILABILITY  OF  ABSTRACT 

0  UNCLASSIFIEOAJNUMITED  □  SAME  AS  RPT.  □  OTIC  USERS 

31.  ABSTRACT  SECURITY  OSSIFICATION  1 

UNCLASSIFIED  | 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Leonard  Abreu 

22b.  TELEPHONE  (Include  Area  Coda) 

(617)  377-4775 

32c.  OFFICE  SYMBOL  1 

GL/OPE  1 

DO  FORM  1473,  M  Mar 


BJ  APR  edition  may  bo  u»d  until  exhausted. 
Ail  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  0*  THIS  PAGE 

UNCLASSIFIED 


i 


smcumTv  cLMwir ic*nm  of  tmi« 


'  approximation  ia  uaed  to  r«duc«  patha  through  the  atmoaphere  to  an  equivalent 

single-layer  path. 


UNCLASSIFIED 

MCUKTY  Ct  >»m>tltaO>  THIS  »AOE 


ii 


TABLE  OP  CONTENTS 


Section 


Page 


1  INTRODUCTION  .  I 

2  THE  MODTRAN  CODE  -  AN  OVERVIEW .  4 

3  MOLECULAR  BAND  MODEL  PARAMETERS  .  7 

3. 1  Line  Center  Parameters  .  7 

3.2  Line  Tail  Parameters .  9 

3.3  Parameter  Data  File . 10 

4  BAND  MODEL  TRANSMITTANCE  FORMULATION  .  13 

4.1  Line  Center  Transmittance .  13 

4.1.1  Curtis-Godson  Approximation  .  15 

4.2  Line  Wing  Absorption .  17 

5  INTEGRATION  INTO  LOWTRAN  7 .  19 

5.1  MODTRAN  SUBROUTINES  .  19 

5.1.1  Subroutine  BMDATA  .  19 

5.1.2  Subroutine  BMOD .  20 

5.1.3  Subroutine  BMLOAD  .  21 

5.1.4  Subroutine  BMTRAN  .  21 

5.1.5  Subroutine  BMERFU  .  21 

5.1.6  Subroutine  BMFLUX  .  21 

5.2  Modifications  to  LOWTRAN  7  .  .....  22 

6  VALIDATION .  2  5 

6.1  Band  Model  Parameters . 25 

6.2  Comparison  with  FASCOD2 . 26 

7  CONCLUSIONS  .  32 

8  REFERENCES . 33 

APPENDIX . 35 


iii 


Accession  Por  • 

~NTIS  GRAAI  “ 
DUC  TAB 
Unannounced 
I  Just  if  1  •ation___ 


|  By - _____ 

!  Distribution/ 

r  -  - - 

!  Availability  Codes 


Dial 


Avail  and/or 
Special 


LIST  OF  FIGURES 


Figure  Page 

1  Atmospheric  Transmittance  for  a  Slant  Path  from  5  to 

10  km  at  15°  from  Zenith  and  Through  the  US  Standard 
Atmosphere  with  No  Haze . .  3 

2  Absorption  from  Lorentzian  Lines  with  Optical  Depths  of 

0.1,  1,  and  10 .  5 

3  Curves  of  Growth  for  Spectral  Lines  with  Combined  Doppler 

and  Lorentz  Line  Broadening .  16 

4  Comparison  of  Continuous  and  Discretized  Triangular 

Slit  Functions . 24 

5  Comparison  of  Band  Model  LOWTRAN  5  and  MODTRAN  Molecular 

Absorption  Coefficients  for  Water  Vapor  at  300  K  .  .  .  .  25 

6  Comparison  of  a  FASC0D2  and  a  MODTRAN  Calculation  ....  26 

7  Comparison  of  a  FASCOD2  and  a  MODTRAN  Calculation 

in  the  Microwave  Region . 27 

8  Comparison  of  Radiances  from  LOWTRAN  7  and  MODTRAN 

( FWHM-20cm~ * ) ,  Both  with  Multiple  Scattering  .  28 

9  LOWTRAN  7  and  MODTRAN  Predictions  for  Radiation  from  the 

H2O  Rotational  Band  for  a  60  km  Limb  Path .  30 

10  MODTRAN  and  SHARC  Predictions  of  H2O  Rotations  for 

60  km  Limb  . . 30 

11  MODTRAN  and  SHARC  Predictions  of  the  15  M®  CO2  Band 

for  60  km  Limb .  3  1 


iv 


1 .  INTRODUCTION 


The  LOWTRAN  computer  codes,  with  LOWTRAN  7^*^  being  th«  biobc  current 
version,  are  widely  used  to  calculate  atmospheric  transmittance  and/or 
radiance  in  the  infrared,  visible  and  near  ultraviolet  spectral  regions; 
LCWTRAN'  7  has  been  extended  to  include  the  microwave  spectral  region.  The 
code  is  easily  used,  runs  quickly  and  provides  the  user  with  a  wide  yatriety 
of  atmospheric  models  and  options.  Its  spectral  resolution  is  20  Full 

Width/Half  Maximum  (FWHM)  with  calculations  being  done  in  5  * 

increments.  _ 

'  This  final  technical  report  describes  ^ork  done  to  increase  LOWTRAN'a 
spectral  resolution  from  20  to  2  cm7* ^  (FWHM).  Specifically,  the  technical 
objectives  for  this  program  were: 

•  to  develop  algorithms  providing  2  cm"1  resolution 
( FWHM ) ; 

•  to  model  molecular  absorption  of  atmospheric 
molecules  as  a  function  of  temperature  and 
pressure  * 

•  to  calculate  band  model  parameters  for  twelve 
LOWTRAN  molecular  species,'  and 

«  to  integrate  the  LOWTRAN  7  capabilities  into  the 
new  algorithms,  maintaining  compatibility  with  the 
multiple  scattering  option.  > 

MGD7T.AN,  the  final  product  of  this  effort,  is  a  moderate  resolution  LOWTRAN 
code.  Molecular  absorption  is  calculated  in  1  cm*'  spectral  bins,  while 
the  otnsi  parts  of  the  calculation  remain  unchanged.  The  molecular  species 
affected  are: 

\  water  vapor,  carbon  dioxide,  ozone,  nitrous  oxide, 
carbon  monoxide,  methane,  oxygen,  nitric  oxide,  sulfur 
dioxide,  nitrogen  dioxide,  ammonia  and  nitric  acid.  ■' 

Tneir  absorption  properties  (the  band  model  parameters)  are  calculated  from 
the  HITRaN  line  atlas, which  contains  all  lines  in  the  0  -  17,900  cm"' 


spectral  region  tint  have  significant  absorption  for  atmospheric  paths. 
The  increased  resolution  of  MODTRAN  spans  the  same  spectral  region; 
calculations  at  larger  wavenumbers,  the  visible  and  ultraviolet  spectral 
regions,  are  performed  at  the  lower  spectral  resolution,  20  cm"-*, 

Creation  of  this  new  code,  MODTRAN,  required  the  development  and 
integration  of  a  completely  new  set  of  band  models  for  calculating 
transmittance.  The  increased  spectral  resolution  is  achieved  using  an 
approach  developed  in  a  previous  effort  which  resulted  in  a  5  cm~* 
option^’^  to  LOWTRAN  5.^5*  In  this  approa-h,  band  model  parameters  were 
calculated  from  the  HITRAN  database  and  used  to  determine  the  equivalent 
width  of  the  absorbing  molecular  gases  in  5  cm-*  spectral  intervals.  In 
the  present  work,  the  band  model  parameters  are  calculated  in  1  cm-* 
intervals.  The  molecular  transmittance  calculation  for  each  bin  has  three 
elements.  First,  the  Voigt  lineshape  o £  an  "average"  line  is  integrated 
over  the  1  cm-*  interval;  when  a  bin  contains  more  than  one  line  of  a  given 
species,  the  lines  are  assumed  to  be  randomly  distributed  with  statistical 
overlap,  finally,  the  contribution  from  lines  whose  centers  are  in  nearby 
bins  is  calculated  as  a  molecular  "continuum"  component,  The  other  LOWTRAN 
components,  which  have  insignificant  spectral  structure  at  1  cm-*,  are 
calculated  at  their  5  cm-*  increments  and  interpolated  to  arrive  at  the 
total  transmittance  for  each  interval.  The  calculation.il  grid  consists  of 
non-overlapping  1  cm-*  bins,  which  are  degraded  to  the  desired  spectral 
resolution  with  an  internal  triangular  slit  function.  Since  these  bins  are 
square  and  non-overlapping,  the  nominal  spectral  resolution  of  MODTRAN  is 
reported  as  2  cm-*  (FWHM). 

The  new  MODTRAN  subroutines  are  written  in  portable  ANSII  standard 
FORTRAN  and  constructed  so  that  their  interfacing  with  LOWTRAN  7  minimizes 
coding  changes.  These  additional  elements  do  not  interfere  with  the 
regular  operation  pf  LOWTRAN  7;  rather  they  represent  an*  additional 
capability  for  higher  spectral  resolution.  All  the  usual  LOWTRAN  options 
like  aerosol  models,  path  options,  multiple  scattering  model, 
user-specified  profiles,  etc.,  have  been  maintained.  The  large  amount  of 
molecular  data  requited  for  the  increased  spectral  resolution  has 
necessitated  using  an  external  data  file  or  band  model  tape. 
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The  comparison  of  MOOT HAN  and  LOWTRAN  7  calculations  shown  in  Figure  1 
illustrates  the  increased  spectral  resolution.  The  figure  shows  the 
transmittance  calculated  at  2  and  20  cm-1  resolution  for  a  low  altitude 
slant:  path  through  the  US  Standard  Atmosphere,  The  2  cm-1  curve  results 
from  the  internal  triangular  slit  function  (IFWHM«2),  and  the  20  cm-1  curve 
is  the  regular  LOWTRAN  7  result  calculated  in  5  cm"1  steps.  The  MODTRAN 
calculation  resolves  line  structure  below  2180  cm-1  (primarily  water),  the 
band  center  of  the  N2O  fundamental  at  2220  cm-1,  and  the  band  center 

at  2284  cm-1. 

An  overview  of  MODTRAN  is  presented  in  Section  2.  In  the  subsequent, 
two  sections,  calculation  of  the  band  model  parameters  and  the 
transmittance  formulation  are  described.  The  changes  to  LOWTRAN  7  and  the 
nev;  MODTRAN  subroutines  are  described  in  Section  5.  In  Section  6, 
illustrative  MODTRAN  calculations  for  initial  validation  are  discussed. 
Conclusions  are  discussed  in  the  final  section. 
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FREQUENCY  (CM-1) 


Figure  1. 


Atmospheric  Transmittance  for  a  Slant  Path  from  5  to  10  km 


at  15°  from  Zenith 

and  Through 

the 

OS  Standard  Atmosphere 

with  No  Haze.  The 

Solid  Curve 

was 

Calculated  with  MODTRAN 

at  2  cm'1  Spectral 

Resolution, 

and 

the  Dotted  Curve  with 

the  Regular  LOWTRAN  7  at  20  cm-1  Resolution. 
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THE  MODTRAN  CODE  -  Mi  OVERVIEW 


2  . 


MODTRAN  is  a  version  of  LOWTRAN  7  with  six  additional  subroutines  that 
provide  the  increased  spectral  resolution.  The  input  data  sequence  for 
hOOTRAN  is  identical  to  LOWTRAN  7's  except  for  two  modifications,  an 
additional  parameter  on  Caids  1  and  A.  A  logical  parameter,  MODTRN,  has 
been  added  to  the  front  end  of  CARD  1 

READ  ( IRD,  ’  (I.1.I4, 12I5,F8.3,F7.  2)  *  )MOnTRN, MODEL,  I  TYPE ,  IEMSCT, 

1  IMULT.Ml ,M2,M3,M4,M5,M6,MDEP, 1M ,  NOPRT,  TBQUND,  SAI.B 

and  the  input  to  CARD  A  has  been  changed  to  integer  format  with  a 
resolution  parameter,  IFWHM,  added 

READ ( IRD, ' (4110) ' ) IV1 , IV2, IDV, IFWHM 

MODTRN  is  simply  a  switch  which  when  set  to  F  (false)  causes  the  tegular 
LOWTRAN  7  to  be  run  and  when  set  to  T  (true)  activates  MODTRAN.  The 
parameter  IFWHM,  which  is  only  read  if  MODTRN  is  true,  specifies  the  full 
width  at  half  maximum,  FWHM,  of  an  internal  triangular  slit  function. 

MODTRAN  and  LOWTRAN  7  differ  in  their  approaches  to  calculating 

molecular  transmittance.  For  several  different  spectral  intervals 
LOWTRAN  7  uses  a  one-parameter  band  model  (absorption  coefficient)  plus 
molecular  density  scaling  functions.  The  MODTRAN  band  model  n-.as  three 
temperature-dependent  parameters,  an  absorption  coefficient,  a  line  density 
parameter  and  an  average  linewidth.  The  spectral  region  is  partitioned 

into  1  cm"  '  bins  for  each  molecule.  Within  each  bin,  contributions  from 
transitions  whose  Tine  centers  fall  within  the  bin  are  modeldd  separately 

from  nearby  lines  centered  outside  of  that  bin.  Figure  2.  The  absorption 

due  to  lines  within  the  bin  is  calculated  by  integrating  over  a  Voigt  line 
shape. The  Curt i s-Godson^ » ® ^  approximation,  which  is  accurate  for  the 
moderate  temperature  variations  found  in  the  earth’s  atmosphere,  is  used  to 
replace  multilayered  paths  by  an  equivalent  homogeneous  one. 
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FREQUENCY  (CM'1) 


Figure  2.  Absorption  from  Lorentzian  Lines  with  Optical  Depths  of 


0.1,  1  , 

and  10.  Halfwidth  was  set  to  0. 

1  cm  1 

.  In 

the 

Band  Model 

Transmittance  Formulation,  Absorption 

from 

Li  ne 

Centers 

(Segment  of  Curves  Falling  Within 

the  Spectra  1 

Bin 

Denoted 

by 

the  Dashed  Lines)  is  Modeled 

Separate  Ly 

from 

Absor pt ion 

Due  to  Line  Tails  (Outside  the 

Dashed 

Li nes ) . 

The  k-distr ibution  method,  which  is  used  in  LOWTRAN  7's  multiple 
scattering  treatment  to  correct  for  averaging  over  large  spectral 
intervals,  is  not  used  in  the  MODTRAN  version.  Because  three 
(monochromatic)  k  values  arc  used  for  the  5  cm”*  steps  of  LOWTRAN  7,  the 
1  cm''1  MODTRAN  steps  provide  on  equivalent  accuracy  for  the  multiple 
scattering  option. 

MODTRAN  is  better  suited  than  LOWTRAN  for  atmospheric  paths  which  lie 
completely  above  30  km.  This  is  due  to  the  integration  over  the  Voigt 
lineshape  combined  with  the  explicit  temperature  and  pressure  dependencies 
of  the  band  model  parameters.  The  Voigt  lineshape  is  necessary  at  these 
altitudes  because  the  Doppler  linewidth  is  greater  than  the  I.orentz.  The 
20  cm"1  versions  of  LOWTRAN  suffer  because  they  use  a  single  set  of  band 
model  parameters  (nominally  sea  level  at  296  K)  coupled  with  spectrally 


is  also 


independent  scaling  functions  for  the  molecular  densities.  It 
worth  mentioning  that,  for  paths  which  lie  completely  above  60  km,  another 
problem  arises:  many  of  the  molecules  are  no  longer  in  local  thermodynamic 
equilibrium  (LTEl) .  This  means  that  che  strengths  of  some  molecular  bands 
can  no  longer  be  determined  from  the  ambient  temperature.  MODTEAN  gives 
reasonable  results  for  those  bands  which  are  in  LTEj  the  problem  is 
identifying  those  spectral  regions  which  are  not  LTE.  For  additional 

information  and  discussion,  the  user  is  directed  to  the  Geophysics 
(9) 


Handbook . 


MOLECULAR  BAND  MODEL  PARAMETERS 


3  . 

The  basic  idea  behind  band  model  techniques^*  *  U  is  to  determine  a 
set  of  parameters  from  which  transmittance  over  finite  frequency  intervals 
can  be  calculated.  In  MODTRAN,  three  hand  model  parameters  are  used,  un 
absorption  coefficient,  a  line  density  and  a  line  width.  The  absorption 
coefficient  measures  the  total  strength  of  lines  in  an  interval.  The  line 
density  is  a  line  strength  weighted  average  for  the  number  of  lines  in  the 
interval,  and  the  line  width  parameter  is  a  line-strength  weighted  average 
line  width. 

MODTRAN  uses  a  bin  width  of  one  wavenumber,  av  »  1  cm'*.  Line  data 
from  the  HITRAN  database, the  AFGL  line  atlas,  is  used  to  calculate  the 
band  model  parameters.  The  tape  contains  data  on  molecular  lines  in  the 


f  l  equem.  y 

range  0  to 

17  900 

cm  * .  For 

each 

molecule  with  lines  whose 

centers  fall  within  a 

given 

spectral 

bi  n , 

the  temperature 

dependent 

absorption 

coet  t  ic  lents 

and 

line  densities 

along  with  the 

line  width 

parameter 

are  stored 

for 

subsequent 

use 

in  calculating 

molecular 

absorption;  a  single  temperature-dependent  absorption  coefficient  parameter 
is  used  to  determine  the  tail  contributions  to  each  spectral  bin  from  lines 
centered  in  nearby  bins. 

In  the  next  two  subsections,  the  calculation  of  the  molecular  band 
model  parameters  is  described,  and  the  formatting  of  the  da  a  file  is 
discussed  in  che  third  subsection.  Formulas  used  to  calculate  the 
parameters  are  given,  along  with  a  discussion  of  their  dependence  on 
temperature  and  pressure. 

3 . 1  Line  Center  Parameters 

Each  frequency  bin  corresponds  to  s  1  cm-1  interval  and  contains 
parameters  for  molecules  with  lines  in  that  interval.  The  molecules  for 
which  band  model  perameters  have  been  determined  are: 

h2q,  co2,  o3,  n2o,  CO,  cha,  o2,  NO,  so2,  NO 2 ,  nh3  &  hno3 
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The  molecular  absorption  coefficient*  (S/d)  (cm”1  amagats”1)  are 

calculated  at  5  reference  temperatures: 

T  -  200,  225,  250,  275  &  300  K  (1) 

Linear  interpolation  is  us  Id  to  calculate  absorption  coefficients  at 
remperatures  between  200  and  300  K.  For  temperatures  below  200  and  above 
300  K,  the  extreme  values,  (S/d) (T»200K)  and  ( S/d ) (T-300K) ,  respectively, 
are  used.  The  absorption  band  model  parameters  are  calculated  from  the 
individual  line  strengths, 

( s/d )  -  1~  y  S j (T)  .  (2) 

j 


Here  S j  (T)  is  the  integrated  line  strength  at  temperature  T  of  the  j'th 
line  of  molecule  m  iri  bin  i.  The  line  strength  at  an  arbitrary  temperature 
is  scaled  from  the  HITRAN  line  strength  at  its  standard  temperature, 
Ts  =  296  K,  by 

qr(Ta)Qv(T8) 
s  j  T  Qr (t)qv (T) 


1  -  exp(-hc v;/kT) 

1  -  exp ( -he v jj kTs ) 


T-Tc 


T  T£ 


Sj(Ts) 


(3) 


where  Qr  and  Qy  are  the  rotational  and  vibrational  partition  lunctions,  Ej 
is  the  energy  of  the  lower  transition  state,  and  vj  ie  the  transition 
frequency.  The  constants  are  the  speed  of  light  (c),  the  Boltzwan  constant 
(k),  and  the  Planck  constant  (h). 

A  collision  broadened  or  Lorentz  line  width  parameter  1°  is  defined  at 
STP  (TQ  ■  273  1  5  K,  PQ  =  1013.25  mbar),  A  single  value  can  be  stored 
because  the  pressure  and  temperature  dependence  of  the  Lorentz  line  width 
is  easily  modeled, 

YC(T.P)  “  --  (T0/T)x  (4) 
*0 
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where  the  exponent  x  has  been  set  to  1/2  for  all  molecules  except  C02 ,  for 
which  x  -  3/4.  The  band  model  parameter  is  calculated  as  a  line 

scrength  weighted  average  over  Che  tabulated  Lorentz  line  widths  "1|'c*j^Ts^ 

-  (Ts/Td)x  [J-rc.jCT*)  Sj(Ts)]  /  r£Sj(Ts)]  .  (5) 

j  J 

Like  the  absorption  coefficients,  the  line  density  band  model 
parameters  (1/d)  (cm)  are  calculated  at  the  five  reference  temperatures  and 
interpolated  when  used  by  the  band  model  subroutines.  The  line  density  is 
defined  by 

N  N 

<*'“>  "  b  <2  sj  >2  '  J  sj2  (6) 

j  =  1  j  “  1 


This  definition  for  the  line  spacing,  which  is  derived  in  the  appendix, 
produces  a  smaller  value  than  the  usual  definition  involving  a  sum  over  the 
square  root  of  the  line  strengths.*'10'11)  The  new  form  results  when 
account  is  taken  of  finite  bin  widths.  The  absorption  of  a  line  within  a 
finite  bin  is  Less  than  its  total  line  strength:  this  is  consistent  with  a 
decreased  value  for  (1/d). 


3 . 2  Line  Tall  Parameters 

The  line  tail  parameters  consist  of  line  contributions  from  lines 
located  outside  of  a  given  bin  but  within  ±  25  cm-1.  The  line  tail 
absorption  coefficient  band  model  parameters  C  (cm-1  amagat”1)  are 
determined  by  integrating  the  Lorentz  line  shape  over  this  interval 


C 


i  +  25 

k=i-25 


+ 


I  (S/d)  7cJk 
(k-i ) 2  +  1/4 

[(S/d)  7c]k 
(k+i)2  +  1/4 


f [ (k~i) Av] 


f[(kti)Av]  ) 


(7) 
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where  che  delta  function  serves  to  exclude  the  k-i  term  from  the  oum  (i.e., 

the  line  center  contribution),  end  f(hvj  is  a  lineahape  for®  factor.  The 

form  factor  is  1.0  within  25  cm"*  of  the  line  centers.  Except  for  HjO  ar.d 
C02,  tail  contributions  beyond  25  cn“*  are  assumed  negligible  and  not 

included.  The  usual  LOWTRAN  7  water  continuum  consists  of  tail 

contributions  from  lines  located  beyond  25  cm”*  plus  extrapolated  (flat) 
values  of  this  contribution  within  25  cm”*  (for  smoothness).  For  CO2 ,  the 
continuum  from  FASC0D2  ^  * ^  >  *^  has  been  added  to  the  C  to  account  for  the. 
tail  contributions  from  lines  beyond  25  cm”*, 

hcv-  Ts 

C  -  C  +  '^tanh  { — — — }  C^)  (8) 

Here,  C(v^)  is  the  frequency  interpolated  value  from  FASC0D2 ‘ s  block  data 
/FC02/.  For  both  H20  and  CO),  the  value  of  C  has  also  been  reduced  by  an 


amount  equal 

CO 

its  value 

at  2  5 

cm”*  from  the  line  center 

since  this 

cont  r ibut ion 

is 

included 

a  1  ready 

in  the  continuum  data. 

The  C  is 

proper t iona 1 

to 

pressure 

(which 

arises  from  the  Lorentz 

line  width, 

Equation  (4)) 

C(P)  - 

P 

*7 

C(PG) 

(9) 

3.3  Parameter  Data  File 


Because  of  the  large  amount  of  data,  the  band  model  parameters  are 
stored  in  an  external  file  that  is  written  in  binary  format;  this  a' lows 
for  quicker  access  during  the  calculation.  Each  entry  corresponds  to  a 
1  cm  *  interval  and  contains  a  molecular  parameter  set.  Since  no  data  for 
molecules  without  lines  in  a  given  interval  are  stored,  (a  parameter 
identifying  the  active  species  is  included. 

The  first  entry  of  a  parameter  set  is  the  bin  number,  i.  From  the  bin 
number,  the  midpoint  of  the  interval  is  calculated 

Vj  »  i  hv  ,  (iO) 
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and  all  lines  whose  centers  fall  in  the  half-opened  interval 
( Vj  -  aw/2,  V|  +  aw/ 2]  contribute  to  bin  i. 

The  molecular  parameter  set  is  identified  by  the  parameter  m.  The 
HITRAN  database*^  convention  is  used  for  this  label 

m  1  23  4  5  6  7  8  91011  12 

molecule  H20  C02  O3  N20  CO  CH4  02  NO  S02  N02  NH3  HNO3 

The.  next  entries  in  the  parameter  set  are  the  molecular  absorption 
coefficients  (S/d)  (era-''  amagats"^)  calculated  at  the  five  reference 
temperatures.  These  entries  are  followed  by  the  STP  Lorentz  half  width, 
nr£,  multiplied  by  10^  and  stored  as  an  integer.  Line  spacing  parameters 
(1/d)  for  the  five  reference  temperatures  complete  the  line  center 
parameter  sets. 

For  line  tails,  each  line  contains  data  on  one  or  two  molecules. 
These  line  tail  parameter  sets  use  the  same  format  as  the  line  center 
parameter  sets.  Again,  the  first  entry  is  the  bin  number  i  and  the  second 
entry  is  the  molecule  designation  m.  To  recognize  that  these  parameter 
sets  denote  line  tail  contributions,  their  molecule  labels  are  offset  by  12 

m  13  14  15  16  17  18  19  20  21  22  23  24 

molecule  H20  C02  03  N20  CO  CH^  02  NO  S02  N02  NH3  HNO3 


The  continuum  parameters,  C,  are  stored  in  place  of  the  (S/d).  Unless  all 
tail  contributions  have  been  defined  for  frequency  bin  i,  the  molecular 


designation  and  continuum  parameters 
on  the  same  parameter  set. 

A  sample  of  the  formatted  data 
shown  below: 


2294 

1 

2.5611-07 

9.8:61-07 

3  1571-06 

8.  694 E  - 06 

2.0831-05 

2294 

2 

1.  3841-rOO 

1  4361+00' 

1  5981+00 

1  9041+00 

2.  3771+00 

2?94 

4 

1.9611-03 

2. 8421-03 

3.  7691-03 

4. 68*1-03 

5.  5401-03 

2254 

13 

1.6811-08 

6  8711-08 

2.0901-07 

5. 1261-07 

1.0711-06 

2294 

15 

5  0521-09 

5.  03/2  -09 

4.8961-09 

4. 6781  09 

4.4161-09 

2294 

18 

3.  3931-06 

3.2941-06 

3  5291-06 

4.0551-06 

4.8161-06 

2295 

2 

3. 1951-02 

7  8651-02 

1.6781-01 

3. 1661-01 

5.  4031-01 

2295 

4 

2  6531-03 

3.8041-03 

5. 0031-01 

6. 1741-03 

7.  2611-03 

2295 

5 

1.5111  06 

l  6031-06 

1.6591-06 

1.6871-06 

1.6911-06 

2295 

13 

2  7961-08 

1  1181-07 

3  4321-07 

8  6281 -07 

1.8601-06 

2295 

16 

3  9031-04 

4  8821-04 

5.  7691 -04 

6  5241-04 

7. 12S1-04 

for  a  second  molecule  follow  the  first 


for  frequency  bins  2294  and  2295  is 


592  1  814E+O0  2.  545E-HJO  3.  273E+0O  3.  3591+00  3  182F+CO 

730  1  I38E+00  1.  '  '5E+O0  1  8371+00  2  4861+00  3.  021E+00 

757  2  2271+00  ■  .oEtOO  2.  5171+00  2. 6271+90  2.  720E+00 

14  8.  3251-02  /  507E  02  7. 245E-02  7.5351-02  8.3561-02 

16  3.  3391-04  4.  363E-04  5  3461  04  0.  230E-04  6.978F-04 

0  0.  000E+O0  0.  OOCE+OO  0.  000£*00  0.  OOOE+OO  0.  OOOE+OO 

688  2.  986E+O0  2  537E+00  2.  3W4E+00  2. 217E+O0  2.  211E+00 

766  2. 729E+00  2  9661+00  3. 16U+00  3.  3201+00  3. 449E+CO 

725  1. 824E+00  l.  898E+00  1.9461*00  1  9751+00  1.991E+00 

14  1  3091  01  1. 2571-01  1.3021-01  1.4441  01  1.675E-01 
18  2.2711-06  3  8341-06  6.  4201-06  9  811E-06  1,3931-05 
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In  bin  2294,  there  are 
data  for  H2O,  C02,  O3, 
data  for  H20,  but  lines 
any  O3  continuum  data. 


H20,  C02  and  H20  line  center  parameter  sets  and  tail 
N20  and  CH4.  In  bin  2295,  there  is  no  line  center 
do  exist  for  CH4.  Also,  bin  2295  does  not  contain 
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BAND  MODEL  TRANSMITTANCE  FORMULATION 


4  . 

4 . 1  Line  Center  Transmittance 

The  band  model  transmittance  f ormu lat ion ^ ^ ^  developed  for  the  5  cm-* 
option  to  LOWTRAM  5^  has  been  used  to  create  a  moderate  resolution  option 
for  LOWTRAN  7.^**  The  expression  used  to  calculate  molecular  transmittance 
is  based  on  a  statistical  model  for  a  finite  number  of  lines  within  a 
spectral  interval,  and  is  given  by 

T  -  d  -  <WsA>/Av)<n>  ,  (11) 

where  t  is  the  transmittance,  <WS  4>  is  the  Voigt  single  line  equivalent 
width  for  the  line  strength  distribution  in  a  spectral  interval,  and  <n>  is 
the  path  averaged  effective  number  of  lines  in  the  bin 

<n>  =  av  <l/d>  .  (12) 

<l/d>  is  the  path  averaged  line  spacing. 

For  large  <n>  [(S/d)  and  av  fixed],  the  transmittance  simplifies  to 

the  more  recognizable  exponential  form.  Beer's  Law,  given  by 

t  -  exp (  -  <WsA>  <l/d>  )  .  (13) 

For  the  relatively  low  temperatures  encountered  in  the  earth's  atmosphere, 
the  number  of  lines  in  a  bin  from  a  single  molecular  species  is  usually 
small  so  that  the  power  law  transmittance  formulation  is  preferred. 

There  are  many  approximations  available  for  calculating  the  equivalent 
width  of  a  Voigt  line  shape;  different  ones  are  valid  for  different 
regimes,  Doppler  or  collision  broadening,  weak  line  or  strong  line,  etc. 
However,  no  single  approximation  is  adequate  for  the  range  of  pressures  and 
optical  path  lengths  encountered  in  atmospheric  transmission  calculations. 
Rather  than  incorporating  different  approximations,  we  directly  evaluate 
the  exact  expression  for  the  equivalent  width  of  a  single  line  in  a  finite 
spectral  interval;  <W sA>  is  given  by 
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m 


<ws4> 


F(X,Y) 


I  -  exp  {-  [Su/d]  Jlnl/rr  F(X,Y)/Ord/d>)dX  ,  (14) 


%S' 

0 

« 

Y  I*  exp  ( -T2 )  dT 

v  J  Y2  +  (X-T)  ^ 


(15) 


xm  -  y  V1"2"  <n>/<-rd/d> 


(16) 


Y  =  ^JTnT  <'Yc/d>/<td/d> 


(17 


\ 

* 


where  F(X,Y)  is  the  Voigt  line  shape  function,  [Su/d]  is  the 
depth,  and  <"Yd/d>  arid  <"irc/d>  are  the  path  averaged  Doppler 
broadened  line  shape  band  model  parameters,  respectively, 
calculate  <Ws4>,  we  separate  its  contributions  as  shown: 


total  optical 
and  collision 
To  accurately 


<w8ft>  -  <ws4°>  -  <ws41> 


<We4°>  ~  ~  fl  -  exp  {-[Su/d]  F(X,Y)/<-Yd/d>)dX 

xm  J 


<w8/> 


Hf 


1  -  exp  {-[Su/d]  VT  nTTif  F(X,Y)/<Td/d>}dX 


m 


(18) 


(19) 


(2d) 


The  tail  contribution,  <Ws4^>  can  easily  be  evaluated  in  terms  of  the  error 
function  since  Xm  >>  Y  for  cases  calculated  with  MODTRAN: 


<Wa41>  o.  exp (- z2 )  +  -Jv  z  erf  (z)  -  1  ,  (21) 

z  ■  — [Su/d ]  <jYc/d)>  /  it  .  (22) 

To  determine  <Ws4®>,  an  interpolation  between  the  Lorentz  and  Doppler 
limits  is  used.  Based  on  an  interpolation  formula  developed  by  Ludwig, 


14 


^■'41  !?>>■■■■"'■  'iHtOTIK  wirwtiNPr 


®t  ,  al.  [Equations  (5-25)  and  (5-26)],^) 
equivalent  widths  are  given  by 


the  Loroncz  and  Doppler 


L  M 


'4  [Su7d]/^Tc/d> 


2  <Td' d>  ,  ...  In2  [Su/d]2, 

D  ■  — *- - rr  ln(  1  +  - - ~) 

ln2  [Su/dJ2  2  <Td/d>2 


(23) 


(24) 


<W,if°>  is  determined  from  Che  following  interpolation  formula  which  is  a 
mote  numerically  stable  form  of  their  formula: 


<Ws«°>2 


2 

lSu/dV  u  -  u-L>  ( i -  d ) /  jr=-rmr= vr mr>  i  .  (25) 


<l/d>*: 


Figure  3  shows  a  comparison  of  Equation  (25)  to  exact  calculations  for  the 
equivalent  width  of  a  single,  isolated  spectral  line  with  a  Voigt 
linesliape.  The  lowest  curve  is  the  pure  Doppler  limit  and  the  highest 
curve  is  the  Lorencz  limit.  The  exact  values  are  taken  from  Penner. 

The  predictions  of  Equation  (25)  are  shown  as  solid  lines  for  the  same 
values  of  the  parameter  Y,  Equation  (17).  The  overall  agreement  between 
Che  two  families  of  curves  illustrates  the  quality  of  the  interpolation 
formula . 


4.1.1  Curtis -Godson  Approximation 

Path  averages  are  calculated  with  the  Curtis-Godson 
approximation. This  approximation  replaces  an  inhomogeneous  path  with 
a  homogeneous  one  by  using  average  values  for  the  various  band  model 
parameters.  The  Curtis-Godson  approximation  is  very  good  for  paths  where 
the  temperature  and  concentration  gradients  are  not.  particularly  steep. 
This  is  the  cade  f6r  atmospheric  paths  where  the  temperature  variations  for 
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( I  n2)  1/2  CSu/d]  /  <'Td/d> 

Figure  3.  Curves  of  Growth  tor  Spectral  Lines  with  Combined  Doppler 
and  Lorentz  Line  Broadening!  The  Exact  Calculations l 
are  Presented  as  Individual  Values.  The  NASA  Formula, 
Equation  (25),  is  Shown  as  Solid  Lines.  The  Values  of  Y 
From  Bottom  to  Top  are:  0,  0.0005,  0.005,  0.05,  0.5,  1.0, 
2,  10,  and  o. 


arbitrary  paths  fall  within  the  range  of  200  tc  300  K.  The  total  optical 
depth  is  a  sum  over  contributions  from  the  individual  layers  and  is  given 
by 

fSu/d]  “  ^  (S/d)4  (Au)s  ,  (26) 
ft 

where  (au)4  is  the'  incremental  absorber  amount  from  layer  A  and  (S/d)A  is 
the  absorption  coefficient  band  model  parameter  at  the  temperature  of  the 
layer  ft.  Note,  frequency  and  species  indices  are  implicitly  assumed  for 
this  and  subsequent  equations. 

The  optical  depth  is  used  as  the  weighting  function  in  calculating  the 
path  averages 
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<l/d> 


(2  7) 


TsOTT  Z 
TsOTT  Z 

ft 

<V-i>  -  fsOTT  Z 

ft 

The  band  model  parameters 
is  the  usual  Doppler  width 

(Va  “  ~  ■JTCKTJWm  ,  (30) 

where  M  is  equal  to  the  molecular  mess  (g/rnole)  of  the  molecule  and  N  is 
Avogadro's  number. 

4 • 2  Line  Wing  Absorption 

The  power  law  transmittance.  Equation  (11),  takes  into  account  only 
lines  which  originate  within  a  given  spectral  interval,  and,  for  these 
lines,  only  that  fraction  of  the  line  profile  which  falls  within  the 
interval  is  included  in  the  computation  of  the  equivalent  width.  This 
approximation  is  reasonable  in  the  strongly  absorbing  region  of  a  band; 
however,  because  the  absorptivity  is  expressed  in  terms  of  the  local  line 
strength  distribution,  it  becomes  a  poor  approximation  in  regions  where  the 
tail  contributions  from  nearby  lines  dominate  the  contributions  from  weak 
or  nonexistent  lines  within  a  given  interval  (bin).  This  typically  occurs 
in  the  center  and  f.ar  wings  of  a  band  (i.e.,  past  the  band  head),  and  also 
in  spectral  intervals  containing  no  lines  which  ore  in  the  vicinity  of 
strong  lines.  For  these  situations,  the  local  absorption  is  dominated  by 
the  accumulated  tails  of  the  stronger  lines  located  outside  the  interval. 
The  effect  of  line  wing  absorption  is  included  in  the  transmittance  by 
adjoining  an  exponential  term 


d/d)  *  (S/d)A  (au)s 


(Vft  n/d)ft  (S/d)ft  (ou)ft  ,  (28) 


(Yd)4  (l/d)A  (S/d ) A  (AU)A  .  (29) 


were  defined  in  the  previous  section,  and  (Tj)^ 
(cm"1) 
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•""•Hnwim  .am 


t  *=  (1  -  <Wuft>/Av)<n>  e“fCul 


(31) 


Where  [Cuj  is  the  total  continuum  optical  depth 


[Cu] 


y  (c)a  (^u) 


A 

The  layer  subscript  ft  on 
temperature  (interpolated) 


ft 


C  labels  both 
dependencies . 


the 


pressure 


(32) 

[Equation  (9)]  and 
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INTEGRATION  INTO  LOWTRAN  7 


Integration  of  che  MODTRAN  subroutines  into  LOWTRAN  7  was  accomplished 
with  minimal  changes  to  the  original  code.  The  interface  between  the 
regular  LOWTRAN  7  and  the  MODTRAN  option  is  made  through  calls  to  two 
subroutines  in  the  LOWTRAN  7  subroutine  TRANS  and  one  call  in  subroutine 
MSRAD.  In  TRANS,  a  single  call  to  subroutine  BMDATA  reads  the  first 
necessary  wavenumber  block  of  band  model  parameters  and  calculates 
wavenumber  independent  quantities.  For  each  wavenumber,  calls  to 
subroutine  BMOD  are  made  once  for  initialization  and  then  additionally  in 
the  loop  over  atmospheric  layers  that  calculates  the  molecular 
transmittance.  In  MSRAD,  the  call  to  FLXADD  is  replaced  by  a  call  to 
BMFLUX  for  the  moderate-resolution  option. 

5 . 1  MODTRAN  Subroutines 

The  MODTRAN  subroutines  are  described  in  this  section.  In  total, 
these  six  subroutines  contain  about  850  lines  of  code. 

5.1.1  Subroutine  BMDATA 

Subroutine  BMDATA  is  called  once  each  calculation  by  subroutine  TRANS. 
This  subroutine  opens  the  binary  band  model  tape,  makes  the  initial  band 
model  tape  reads,  and  calculates  wavenumber -independent  quantities  for 
subsequent  use  by  subroutine  BMOD. 

After  opening  the  band  model  tape,  BMDATA  reads  the  header 
information,  advances  the  tape  to  the  wavenumber  block  containing  the 
initial  frequency,  l'Vl,  and  reads  that  block  of  data  into  common  /BMDCOM/. 
If  the  requested  spectral  interval  is  totally  outside  of  the  band  model 
tape  range,  the  program  simply  performs  a  standard  LOWTRAN  7  calculation. 

For  each  atmospheric  layer,  wavenumber  independent  quantities  are 
calculated  in  BMDATA  and  stored  in  the  appropriate  arrays  for  later  use  in 
BMOD.  The  quantities  stored  are  VT7T7  in  arrays  WT  and  WTS,  and 
temperature  interpolation  indices  in  arrays  JJ,  FF,  JJS  and  FFS.  Arrays 
affixed  with  the  letter  'S'  define  parameters  for  the  solar  paths. 
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5.1.2  Subroutine  BMOD 


Subroutine  BMOD  is  the  primary  focus  of  the  moderate  resolution 
subroutines.  In  BMOD,  the  Cur t is-Godson  sums  are  calculated  (Subsection 
4.1.1),  subroutine  BMTRAN  is  called,  and  the  continuum  contributions  due  to 
tails  of  lines  originating  outside  the  1  cm-^  interval  are  determined.  An 
initialization  call  is  made  every  wavenumber  from  subroutine  TRANS,  and 
subsequent  calls  are  made  for  each  layer  to  calculate  the  transmittance. 
On  the  first  call,  BMOD  retrieves  band  model  data  through  calls  to  BMLOAD, 
initializes  the  Curtis-Godson  sums,  determines  the  frequency  dependent 
Doppler  line  widths.  Equation  (30),  and  reads  the  next  frequency  block  of 
band  model  parameters  if  necessary  (IP  is  the  counting  variable),  On 
subsequent  calls,  the  Curt is-Codson  approximation  is  used  to  calculate 
optical  depth  (SDSUM) ,  line  density  (ODSUM),  Lorentz  half  width  (COLSUM) 
and  Doppler  half  width  (DOPSUM).  The  contributions  from  tails  of  lines  are 
also  calculated  within  the  layer  loop.  Since  these  toils  have  a  smooth 
spectral  structure,  they  form  a  continuum  component,  and  the  simple 
exponential  form.  Equation  (31),  is  used  to  calculate  the  transmittance. 
These  tail  contributions  and  the  molecular  components  from  BMTRAN  are 
combined  (multiplied  together)  and  stored  in  the  matrix  TX  for  subsequent 
use  by  subroutine  TRANS. 

The  calls  made  by  TRANS  to  BMOD  for  calculating  the  transmittance  is 
embedded  in  a  loop  over  atmospheric  layers.  When  only  the  transmittance  is 
calculated,  the  DO-loop  in  TRANS  consists  of  one  pass.  When  atmospheric 
radiation  is  calculated,  the  TRANS  loop  is  over  all  layers,  be  ause  the 
radiation  depends  on  the  incremental  change  in  transmittance  for  each 
layer.  The  Curtis  Godson  approximation  for  the  statistical  band  model 
requites  the  suta  over  layers  for  transmittance  calculations.  In  keeping 
with  the  philosophy  of  minimizing  changes  to  the  basic  LOWTRAN  7  program, 
this  layer  loop  is  done  in  BMOD  for  transmittance-only  calculations  and  in 
TRANS  for  radiation  calculations.  Thus,  the  ’  ver  loop  in  BMOD  consists  of 
just  one  pass  when  TRANS  is  looping  over  all  layers  and  vice-versa. 
Subroutine  BMOD  determines  whether  the  loop  in  TRANS  is  over  one  layer  or 
all  layers  and  then  adjusts  its  own  loop  variable  accordingly. 
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5.1.3  Subroutine  BMLOAD 

Subroutine  BMLOAD  (called  by  BMOD)  loads  band  model  data  for  a  single 
parameter  set  into  the  matrices  SD,  OD  and  ALFO. 

5.1.  4  Subroutine  BMTP.AN 

BMTRAN  is  called  by  subroutine  BMOD  to  calculate  finite  bin  molecular 
t ransmit tances .  The  curve-of-growth  used  in  the  statistical  band  model  is 
based  on  the  equivalent  width  of  a  single  average  line  in  a  1  cm"* 
interval.  Band  model  parameters  are  calculated  in  BMOD  for  an  equivalent 
homogeneous  path,  using  the  Curtis-Godson  approximation. 

As  discussed  in  Subsection  4.1,  the  Voigt  line  shape  is  integrated 
over  the  interval.  When  the  optical  depth  is  less  than  .001,  the  weak,  line 
limit  is  used.  Otherwise,  the  Voigt  line  shape  is  calculated  from 
Equations  (18)  -  (25). 

5.1.5  Subroutine  BMERFU 

BMERFU  is  called  by  BMTRAN  to  calculate  the  finite  bin  width 
correction  to  the  Voigt  integral,  Equations  (21)  and  (22).  For  small 
non-negative  z,  the  error  function  is  calculated  from  the  rational 
approximation  given  by  Abramowitz  and  Stegun,  Equation  (7.1.26). For 
large  z,  the  error  function  is  determined  as  a  continued  fraction  derived 
from  the  asymptotic  expansion,  Equation  (7.1.23)  of  Reference  16. 

5.1.6  Subroutine  BMFLUX 

BMFLUX  is  the  moderate  resolution  version  of  the  LOWTRAN  >7  subroutine 
FLXADD,  Called  by  MSRAD,  BMFLUX  is  vied  to  sum  optical  thicknesses,  Co 
compute  the  diffuse  thermal  and  solar  flux  contributions  from  a  two-stream 
approximation,  and  to  combine  these  fluxes  in  an  adding  routine  which 
determines  total  upward  and  downward  fluxes  for  each  layer.  The 
k-distribution  method  is  not  used  in  this  fast  version  of  FLXADD.  Since 
the  lectial  intervals  in  the  MODTRAN  calculations  are  only  1  ctn“*,  it 
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treats  these  transraittances  as  monochromatic  for  the  multiple-scattering 
fluxes.  Thus,  the  molecular  optical  depth  of  a  layer  is  calculated  in 
MSRAD  as  a  logarithm  of  molecular  transmittances. 

5 . 2  Modifications  to  LOWTRAN  7 

Modifications  to  LOWTRAN  7  have  been  kept  to  a  minimum.  As  mentioned 
in  Section  2,  the  switch  MODTRN  has  been  added  to  /CARD1/.  Only  if  MODTRN 
is  .TRUE,  are  any  of  these  changes  activated. 

Most  of  the  routines  from  LOWTRAN  7  remain  unchanged.  A  number  of 
routines  have  been  modified  only  in  that  the  blank  common  along  with  the 
labeled  commons  /CARD1/,  /CARD4/,  /SOLS/  and  /TRAN/  have  been  changed. 
These  routines  are  ABCDTA,  AEREXT,  AERNSM,  CIRR18,  CIRRUS,  CLDPRF,  DESATT, 
EQULWC,  EXABIN,  FLXADD,  LAYVSA,  PHASEF,  RDEXA,  RDNSM,  RFPATH,  SSRAD  and 
VSANSM . 

A  number  of  routines  have  undergone  minor  modifications.  Routines  GEO 
and  SSGEO  were  altered  in  order  to  make  the  following  changes/additions: 

•  The  matrices  WPATH  (WPATHS)  and  WPMS  (WPMSS)  have 
been  combined  into  the  single  matrix  WPATH 
(WPATHS) , 

•  Layer  pressures  and  temperatures  are  stored  for  use 
by  the  band  model  routines,  and 

•  Curt is-Godson  averaged  pressures  and  temperatures 
are  determined  for  solar  paths. 

STDMDL  calculates  actual  rather  than  scaled  molecular  column  densities  when 
MODTRN  is  .TRUE..  Finally,  routine  MSRAD  computes  molecular  optical 
thicknesses  and  calls  routine  BMFLUX. 

Significant  changes  were  made  to  the  MAIN.  It  has  been  split  into  two 
routines.  The  new  MAIN  consists  of  almost  1000  lines  of  introductory 
comments  and  a  single  call  to  the  new  subroutine  DRIVER.  DRIVER  is  the 
driver  for  MODTRAN  and  it  contains  all  the  executable  statements  from 
LOWTRAN  7's  MAIN.  In  addition,  it  defines  a  pointer  array  called  KPOINT 
that  maps  the  HITRAN  molecular  labels  (1-12)  into  Che  LOWTRAN  7  labels. 
Also,  DRIVER  checks  the  spectral  inputs.  For  moderate  resolution,  the 
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variables  IV1,  1V2  and  1DV  need  not  be  multiples  of  5,  and  the  additional 
variable  IFWHM  must  be  read. 

Considerable  modifications  were  required  for  TRANS.  When  MODTRN  is 
.TRUE.,  TRANS 

•  sets  the  internal  frequency  step  size  to  1  cm-*, 

•  calls  the  MODTRAN  subroutines  to  calculate  molecular 
transmittance , 

•  interpolates  transmitt Ances  calculated  at  5  cm-*  intervals 
for  the  aerosols  and  molecular  continua,  and 

•  employs  a  discretized  triangular  slit  function  with  FWHM  set 
to  IFWHM  to  automatically  degrade  the  1  cm-*  bin  results  to 
the  requested  spectral  resolution. 

In  Subsections  5.1.1  and  5.1.2,  the  calls  to  the  MODTRAN  subroutines  were 
discussed.  Note  that  the  discretized  triangular  slit  function  when  IFWHM 
is  set  to  1  cm-1  is  a  1  cm-*  rectangular  slit  function,  as  shown  in  Figure 
A.  The  figure  also  demonstrates  that,  when  the  printing  step  size 
parameter  IDV  (=2  here)  is  larger  than  IFWHM,  the  frequency  range  is  under 
sampled . 
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Figure  4,  Comparison  of  Continuous  and  Discretized  Triangular  Slit 
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6.  VALIDATION 


Tests  were  performed  to  ensure  the  band  model  parameters  were 
generated  correctly,  and  MODTRAN  transrai ttances  and  radiances  were 
validated  by  comparing  to  FASC0D2  line-by-line  calculations.  A  few 
comparisons  follow. 

6 . 1  Band  Model  Parameters 

The  1  cm”^  band  model  parameters  were  checked  against  the  moderate 
resolution  LOWTRAN  5^,4)  g  cni~*  values.  A  straightforward  degradation  of 
the  newer  higher  resolution  parameters  was  performed.  In  Figure  5, 
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Figure  5.  Comparison  of  Band  Model  l.OWTRAN  5  (Dotted  Curve)  and 
MODTRAN  (Solid  Curve)  Molecular  Absorption  Coefficients 
for  Water  Vapor  at  300  K. 


absorption  coefficients  from  both  sets  of  parameters  before  degradation  are 
plotted.  When  degraded,  the  1  cm*1  parameters  are  indistinguishable  from 
the  curve  with  the  LOWTRAN  5  band  model  parameters. 

6 . 2  Comparison  with  FASCOD2 

In  Figure  6,  MODTRAN  and  FASC0D2  transmittances  are  compared  for  a  low 
altitude  slant  path  through  the  US  Standard  Atmosphere  in  the  spectral 
region  2160  -  2260  cm"1.  The  curves  were  generated  with  a  2  cm-1 
triangular  slit  function  (FWHM).  The  absorption  in  this  region  is 
primarily  due  to  the  N20  band  centered  at  2224  cm"1  and  the  P  branch  of  the 

4.3  CO2  band.  Although  the  agreement  is  very  good,  discrepancies  do 
exist.  The  dip  at  2210  cm"1  stands  out.  The  N20  lines  in  this  region  are 
spaced  about  0.8  to  0.9  cm"1  apart,  and  therefore  most  1  cm"1  bins  contain 
just  one  line.  The  bin  at  2210  cm"1  is  an  exception  containing  two  lines. 
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Figure  6.  Comparison  of  a  FASC0D2  and  a  MODTRAN  Calculation.  N?0 
and  COy  are  the  Primary  Absorbers. 
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resulting  in  che  dip.  Similar  analysis  can  be  applied  to  other 

discrepancies  and  implies  that  better  agreement  results  with  some  further 

degradation  of  the  spectral  resolution. 

Since  che  band  model  parameters  were  generated  down  to  0  cm~*s  MODTRAN 

can  be  used  to  perform  molecular  absorption  calculations  in  the  microwave 

region.  Results  for  a  0.1  km  sea  level  path  through  the  US  Standard 

Atmosphere  are  compared  to  FASC0D2  in  Figure  7.  The  dominant  absorber 

throughout  this  region  is  water.  The  calculations  compare  well. 

As  discussed  in  Section  2,  the  LOWTRAN  7  and  MODTRAN  implementations 

of  the  multiple  scattering  routines  are  different.  LOWTRAN  uses  the 

k-distr ibution  method  to  approximate  the  multiple  scattering  contribution 
_  1 

to  each  5  cm  1  interval  as  the  sum  of  three  monochromatic  calculations. 
The  basic  idea  is  to  partition  the  range  of  absorption  coefficients 
occurring  in  the  5  cm-1  interval  into  three  representative  values,  each 
weighted  according  to  the  size  of  the  sub-interval  it  represents.  Since 
the  intervals  in  MODTRAN  are  only  1  cm-1  wide,  the  partitioning  into 


ca 


FREQUENCY  (CM’1) 


Figure  7.  Comparison  of  a  FASC0D2  and  a  MODTRAN  Calculation  in  the 
Microwave  Region.  The  Calculation  was  Performed  Using  a 
0.1  km  Horizontal  Path  at  0  km  Altitude  and  with  No  Haze. 


27 


sub-intervals  has  been  eliminated,  and  a  single  average  absorption 
coefficient  is  used  for  each  bin. 

Figure  8  shows  a  comparison  of  a  LOWTRAN  7  multiple  scattering 
calculation  (dotted  line}  to  a  MODTRAN  multiple  scattering  calculation 
(solid  line)  from  .6  to  .7  pm,  each  at  20  cm-*  resolution  and  with  a  step 
size  of  5  cm-1.  The  atmospheric  profile  is  the  U.  S.  Standard,  and  the 
calculation  is  performed  using  the  rural  extinction  atmospheric  haze  model 
with  a  ground  visibility  of  5  km.  The  observer  is  at  an  altitude  of  20  km 
looking  straight  up,  and  the  sun  is  60°  from  zenith.  For  this  calculation, 
thermal  radiance  is  negligible  and  the  contribution  of  single  scattering 
(dashed  line)  is  only  about  15%  of  the  total  signal;  the  radiance  is 
dominated  by  multiple  scattering.  Agreement  between  the  two  calculations 
is  excellent,  the  main  difference  being  the  increased  spectral  structure  in 
the  MODTKAN  calculation. 


FREQUENCY  (CM1) 

Figure  8.  Comparison  of  Radiances  from  LOWTRAN  7  and  MODTKAN 
( FWHMa20cm~  1 ) ,  Both  with  Multiple  Scattering.  Visibility 
is  5  km  and  t~he  Sun  is  60^  from  Zenith.  The  Observer  i s  at 
20  km  Altitude  Looking  Straight  up 
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LQWTRAN  ?  and  MODTRAN  also  differ  in  their  approaches  to  handling 
molecular  transmictance.  Since  the  LOWTRAN  model  has  been  optimized  for 
296  K,  sea  level  paths,  LOWTRAN  should  not  be  used  for  atmospheric  paths 
completely  above  30  km.  MODTRAN,  on  the  other  hand,  uses  a  Voigt 

lineshape,  which  is  applicable  at  higher  altitudes.  Figures  9-11 
demonstrate  MODTRAN's  high  altitude  capabi  1  i t ies .  First,  Figure  9  shows 

that  LOWTRAN  7  and  MODTRAN  do  indeed  predict  vastly  different  radiances  at 
higher  altitudes.  Radiation  levels  from  H2O  rotations  along  a  60  km  limb 
path  are  shown.  The  LOWTRAN  spectral  radiances  are  much  too  low  at  these 
altitudes.  To  demonstrate  that  MODTRAN  calculations  are  correct, 

validations  have  been  performed  against  SHARC,  the  Strategic  High 

Altitude  Radiation  Code.  SHARC  perforins  NLTE  (non-local  thermodynamic 
equilibrium)  calculations  from  60  to  300  km  altitude.  However,  at  60  km, 
vibration  state  populations  are  essentially  LTE  and  H2O  rotations  are 

always  treated  as  LTE  in  SHARC,  so  comparisons  between  MODTRAN  and  SHARC 
should  produce  similar  results.  With  a  60  km  limb  path,  the  two  codes 

predict  similar  spectral  radiances  for  H2O  rotations  (Figure  10)  and  for 
the  15  jjm  CO2  band  (Figure  11),  which  is  mosrly  LTE.  The  SHARC 
calculations  were  done  at  a  spectral  resolution  of  0.5  cm-1  and  degraded  to 
1  cm-1  (FWHM). 
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LOWTRAN  7  and  MODTRAN  Predictions  for  Radiation  from  the 
HgO  Rotacional  Band  for  a  60  km  Limb  Path.  Ac  These 
Altitudes,  LOWTRAN  Under  predicts  the  Radiance. 


FREQUENCY  (CM"1) 

Figure  10.  MODTRAN  and  SHARC  Predictions  of  H2O  Rotations  for  60  km 
Limb. 
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7.  CONCLUSIONS 


A  radiative  transfer  model  with  a  spectral  resolution  of  2  cm'*  (FWHM) 
has  been  developed  based  on  the  LOWTRAN  7  model.  The  new  code  has  been 
given  the  name  MODTRAN  (Moderate  resolution  LOWTRAN).  Initial  validation 
studies,  based  on  comparisons  to  FASC0D2,  indicate  no  significant 
discrepancies  at  2  cra~*  resolution,  FWHM.  MODTRAN  can  be  used  to  calculate 
atmospheric  radiance  at  upper  altitudes,  specifically  for  any  path  in  which 
the  LTG  (local  thermodynamic  equilibrium)  approximation  is  valid. 
Molecular  absorption  is  modeled  as  a  function  of  both  temperature  and 
pressure  for  twelve  molecular  species  -  water  vapor,  carbon  dioxide,  ozone, 
nitrous  oxide,  carbon  monoxide,  methane,  oxygen,  nitric  oxide,  sulfur 
dioxide,  nitrogen  dioxide,  ammonia  and  nitric  acid.  MODTRAN  also  has  a 
multiple  scattering  option  which  enables  the  code  to  calculate  both  thermal 
and  solar  multiple  scattering  contributions.  A  final  version  of  the  code 
has  been  delivered  to  the  Geophysics  Laboratory. 
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APPENDIX 


In  MODTRAN,  che  molecular  transmittance,  t,  from  lines  originating 
within  a  spectral  bin,  Av,  is  determined  from  an  expression  of  the  form 

Av/2 


T  ■  (  e" 

Sub ( v) 

dv  )n 

(A-l) 

0 

where  b(v)  is  a  line 

shape 

function,  u  is  che 

absorber 

amount,  and  S 

and  n 

are  functions  of  the 

absorption  coefficient 

(S/d)  and 

1 ine  density 

(1/d) 

band  model  parameters 

s  -  <s/d! 

(l/d) 

» 

(A-2) 

n  “  (1/d)  Av 

(A-3) 

Equation  (A-l)  compare  to  Equations  (11)  and  (14). 

A  standard  method  for  determining  band  model  parameters  was  laid  out 
by  Goody.  jpe  weak  and  strong  line  limits  of  a  two  parameter 

transmittance  expression  are  determined.  These  results  are  equated  to  the 
weak  and  scrong  line  limits  of  che  Ladenburg-Reiche  function,  and  che 
resulting  equations  are  solved  for  the  unknown  parameters.  When  Goody 
applied  his  approach  to  two  trial  transmittance  formulas,  both  gave 
expressions  for  the  band  model  parameters  of  the  form 

U-4) 

N 

(l/d)  *  ±-  <£  V5p2  /  7  si  (A-5) 

i»  1 


(S/d) 


-  y 

Z* 

i?i 
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Here,  N  is  che  number  of  lines  and  is  the  integrated  line  strength  of 
line  i  (the  line  spacing  formula  assumes  all  N  lines  have  the  same  half 
width) . 

It  is  worth  noting  char.  Goody's  approach  is  not  quite  appropriate  for 
the  transmittance  expression  in  MODTRAN.  Equation  (A-l)  models  the 
transmittance  through  a  set  of  lines  whose  centers  fall  within  a  given 
finite  spectral  bin.  The  Ladenburg-Reiche  function,  on  the  other  hand, 
models  the  total  absorption  due  to  those  lines  -  not  just  their 
contribution  within  a  finite  spectral  bin. 

Instead  of  equating  the  weak  and  strong  line  limit  of  Equation  (A-l) 
to  the  Ladenburg-Reiche  function,  we  equate  the.  weak  and  strong  line  limit 
of  Equation  (A-l)  to  Plass'  expression ^ 2^)  for  the  transmission  due  to 
lines  whose  centers  are  randomly  distributed  in  a  given  spectral  interval, 


Av/2 


Av/2 


Av 


J 


exp[-Sub(v)]  dv  }n  2k  { 


Av 


/ 


exp[-S£Ub(v) ]  dv  )  , (A-6) 


Here,  all  N  lines  are  assumed  to  have  the  same  half  width,  but  the  line 
shape  function  is  arbitrary. 

If  Bj  is  defined  by  the  equation 


Av/2 


0 


dv 


(A-  7 ) 


then,  in  the  weak  line  limit,  Equation  (A-7)  becomes 


1  +  u  (Bj  n  S),+  u  2 


(  n  S  )2 


n  S2  ]  +  0(u3)  - 


N 


+ 


u 


2 


+  0(u3)  .(A-8) 


The  first  order  equation  in  u  reestablishes  Equation  (A-4) 
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n  S 


(A- 9) 


To  determine  the  strong  line  limit  of  Equation  (A-6),  we  assume  b(v) 
to  be  a  non-increasing  bounded  positive  function  on  [0,«).  Consider  a 
general  integral  of  the  form 


I  “ 


Av/2 


2 

av 


/ 


,-Sub ( v)  dv 


(A- 10) 


It  follows  that 


0  <  e-Sub(0)  <  i  <  t>-Sub  (av/2  ) 


(A- 11) 


This  inequality  requires  that  I  falls  off  exponentially  with  large  u  and 
that  the  exponent  is  linear  in  u.  The  exponential  decay  factor,  C,  is  given 
by  the  limit  as  u  -*  «  of  the  negative  logarithm  of  I  divided  by  u 

C  -  -  lira  [ln(I)/u]  .  (A- 12) 

Applying  L'hopitai's  rule,  one  obtains 


C  »  S  <b(v)> 


(A-13) 


where 


<b(v)> 


lim 


Av/2  Ckv/'i 

b ( v)  e-Sub(v)  d v  j  ^  e~ Sub(v)  dv 


f 


(A-14) 


Note  that  this  limit  is  independent  of  the  value  of  S  as  long  as  S  >  0. 

Taking  the  logarithm  of  both  sides  of  Equation  (A-6),  dividing  by  u, 
and  equating  the  asymptotic  limits,  one  obtains 

N 

n  S  <b(v)>  ■  Y  S:  <b(v)>  .  (A-15) 

i-1 
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Since  <b(v)>  is  line  independent,  Che  surprising  result  is  chat  Equation 
(A-9),  che  weak  line  limit,  is  also  the  strong  line  limit! 

To  determine  n  and  S,  we  return  to  the  weak  line  limit,  Equation 
(A-8),  and  consider  the  equation  obtained  from  the  coefficient  or'  the 
quadratic  terms  in  u.  This  equation  gives 


(A-16) 


The  resulting  expressions  for  the  band  model  parameters  are 


(S/d) 


(A-4) 


(1/d) 


(A-17) 


(S/d)  has  the  same  form  as  before,  but  a  smaller  value  for  (1/d)  is 
obtained . 


